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c5 1 Introduction 

The discovery of Hawking radiation0 has not only solved some contradictions inside 

black hole thermodynamics but also revealed a profound intrinsic relations among quantum 

mechanics, thermodynamics and general relativity. It is an important subject of black 

hole physics to investigate the thermal properties of various black holes, o Thanks to the 

existence of various astrophysical processes such as evaporation and accretion, @ a black 

hole in the universe must change with time, thus to make a deeper investigation on the 

Hawking radiation of an evaporating black hole is of much significance to understand 

such dark celestial bodies as black holes. In recent years, researches on the Hawking 

radiation of non-stationary black holes attract much more attention, but most of these 

studies concentrated on the quantum thermal effect of scalar fields. The investigation 
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of the Hawking radiation of fermions in the non-static black holes is merely limited to 
those in spherically symmetric ones.0 Recently, by making use of the generalized tortoise 
transformation proposed by Zhao Zheng, Wu and Caii has successfully dealt with 
the Hawking radiation of Dirac particles in a non- stationary Kerr black hole.S They 
simultaneously treated the first-order and second-order Dirac equations and derived both 
the event horizon equation and the Hawking temperature which coincide with previous 
results. A New quantum effect observed in the thermal radiation spectrum of fermions 
probably originates from the coupling of the spin of particles with the angular momentum 
of black holes. It is absent from the Bosonic radiation spectrum of scalar particles. @ 

In this paper, we use the method developed in Ref. ||] to cope with the Hawking 
radiation of massless fermions namely Weyl neutrinos in a rectilinearly non-uniformly 
accelerating Kinnersley black hole.0 The corresponding quantum thermal effect of scalar 
particles in the same space-time has already been studied in Ref. ||. 

2 The Weyl equation in the background space-time 

The metric of a rectilinearly non-uniformly accelerating Kinnersley black hole that we 
will take into account can be rewritten in the advanced Eddington-Finkelstein coordinate 
system as'0 

ds 2 = 2dv(Gdv -dr- r 2 fd6) - r 2 [dQ 2 + sin 2 Od<p 2 ) , (1) 

where 2G = 1 — 2M/r — 2arcos6 — r 2 f 2 , f = —asm 6, the parameter a = a(v) is the 
magnitude of acceleration, the mass M[v) of the hole is a function of time v. 

We choose such a complex null-tetrad {l,n, m, m} that satisfies the orthogonal con- 
ditions In — —m • m = 1. Thus the covariant one- forms can be written as 

I = dv , n = Gdv — dr — r 2 fd6 , 
— r 

m = —j= (d9 + i sin 9dip) , 
v2 

— T 

m = —= (d9 — i sin 6 dip) , (2) 
and their corresponding directional derivatives are 



D = —d r , A = d v + Gd r , 
V2r V smv ■ 1 

-F 1 / Q % 



\[2r V sin 6* 



r'fd r + d e - —d v . (3) 
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It is not difficult to compute the non- vanishing Newman-Penrose (NP)Q spin coefficients 
as follows (we denote G >r = dG/dr) 

p=l/r, fx = G/r, 7 = -G r /2, r = -jt = f/V2, /3 = cot 0/(2 V%r) , 

a = -cot6/(2V2r) + f/y/2, u = [(2rG - r 2 G,.)f + r 2 f tV + G >e ] /(V2r) . (4) 

Inserting the related spin coefficients and directional derivatives into the spinor form 
of the Weyl neutrino equation's in the NP formalism 

(D + e- p)T)i -(S + 7T- a)r] = , 

(A + /i- 7 >7o-($ + /3-T)77i = 0, (5) 



we obtain 



V lVl + -—(£- r 2 fV 2 ) 7/o = 



V2r 

(d v + GV X + G >r /2) m - (£ f - r 2 fV,) m = 0, (6) 
in which we have defined operators 



1 % 1 % 
V n = d r + n/r , C = d e + - cot 9 r~xd v , C ] = d e + - cot 9 + -— -d v 

2 sin 9 2 sm9 ' 



With further substitutions Xi = V2 r Vii Xo = Vo? Eq. ©, we can reduce them to 
the form 



V oX i + (C - r 2 fV 2 ) Xo = 0, 

r 2 {2d v + 2GVi + G, r ) xo ~ (C ] - r 2 fV ) X i = . (7) 



3 Event horizon equation 

Though Eq. (^) can not be decoupled with respect to 9, however to deal with the 
problem of Hawking radiation, one should be concerned about its asymptotic behaviors 
near the horizon only. Basing upon of the symmetry of the space-time under consideration, 
we introduce a generalized tortoise coordinate transformation as follows 

1 

r* = r + — ln[r - r h (v,9)] , = v - v , 9* = 9 - 9 , (8) 
Zk 



namely, 



dr* = dr H ; (dr — r h dv — r' h dff) , dv* = dv , d9* = d9 . 

2K[r-r h ) 
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where = Th(v,9) is the location of the event horizon, k is an adjustable parameter 
to be determined, the parameters vq and 8q are arbitrary constants, all of them remain 
unchanged under the tortoise coordinate transformation. r h = d v r h and r' h = dgr h can be 
viewed as parameters depicting the evolution of the event horizon. All induced relations 
among the first-order and second-order derivatives from the transformation (El) can be 
referred to Ref. 0] and are not listed here. 

Now we apply the transformation (EJ) to Eq. (|7|) and take the r — > r h (vQ,9 ), v — * v 
and 8 —>■ 6q limits, then we can obtain its limiting form near the event horizon as 

(< + rlf) £*i + 2r\ (G - f h ) = 0. (9) 

Obviously the condition that Eq. (|9]) has a nontrial solution for xo an d Xi i s that its 
determinant vanishes, which gives the following equation to determine the location of the 
event horizon 

2G-2r h +(r h f+ T M 2 = 0. (10) 
v r h > 

It is easily to see that this result is consistent with that inferred from the null-surface 
condition: g^diFdjF = 0. Making a similar prescription as done in the above to the 
null-surface equation 

(2G + r 2 f){d r Ff + 2d v Fd r F - 2fd e Fd r F + \(d 6 F) 2 + -^-(d^F) 2 = , (11) 
we can arrive at an equation 




2G + r 2 h f 2 -2r h + 2fr' h + ^ T — F =0, 



which results in the same one (|10|) that determines the event horizon after we let the 
coefficient in the front of the derivative ^^-FJ to be zero. Eq. (fl0|) shows that 
depends not only on v but also on 9. This means that the location of the event horizon 
and the shape of the black hole change with time. 



4 Hawking temperature 

To investigate the Hawking radiation of Weyl neutrinos, we are now in a position 
to consider the second-order form of the Weyl equation. After a direct calculation to the 
second-order Weyl equations 

[r 2 V 2 (2d v + 2GV X + G, P ) + {O - r 2 fV 2 )(C - r 2 fV 2 )] Xo = , (12) 
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and 



[r 2 {2d v + 2GV X + G r )X? + (C - r 2 fV )(tf - r 2 fV )] Xl 
= 2r 2 {[{2rG - r 2 G, r )f + r 2 f, v + G, e ]d r 
+ (rG, r -3G + r 2 G rr /2) / + 2rf\ v + G fi /r + G >r e/2} Xo 



(13) 



we can get their explicit forms as follows 

\2r 2 G + + 2r 2 d 2 vr + d 2 + -\-d 2 - 2r 2 fd 2 e + 

1 sin v * sin v 

+ (cot 9 - Arf)d e + 4rd v + (3r 2 G r + 6rG + 6r 3 f 2 - 2r 2 f cot 0)0, 



+r 2 G rr + ArG r + 2G + 6r 2 f 2 - Arf cot 9 



1 



1 



4 sin 2 9 4 



Xo = 0. 



(14) 



and 



(2r 2 G + r 4 f 2 )d; + 2r 2 d 2 r + 4 . 

sin 9 r 

+ cot + (r 2 G r + 2rG + 2r 3 f 2 - 2r 2 / cot 0)9, 
= 2r 2 {[(2rG - r 2 G, r )f + r 2 f\ v + G fi ]d r 



l —d 2 v - 2r 2 fd 2 , - 

sin 9 
1 1 



4 sin 2 9 4 



Xi 



+ (rG, r -3G + r 2 G rr /2) / + 2rf, v + G fi jr + G, re /2} Xo • 
Application of a similar procedure as in the last section to Eqs. (]14,15) leads to 



(15) 



^ + 2r 2 (2G-r h ) + 2r^ 2 + 2/r 2 r; 



|^Xo - 2 (jr\ + r' h 



& 2 



dr*d6. 



-Xo 



+2ri 



d 2 



dr^dv* 



Xo + 



- A + 3r 2 G r + r h (QG - Ar h ) + Qr 3 J 2 



d 



-2r 2 J cot 9 + (4fr h - cot 9 )r' h - r£ — Xo = 



and 



A + 2r 2 (2G - r fc ) + 2rJ/ a + 2fr 2 h r' h 



dr 
9 2 



(16) 



arf 



rr + r, 



^2 



Xi 



a 2 



dr*dv* 
d 



-A + r\G, r + 2r h G + 2r\f - 2r\f cot 9 - cot 9 r' h 



Xi = 2r 2 h (2r h G-r 2 h G, r )f + r 2 h l v + G, e 



dr 

r 'h + rlf 
G-r h 



_d_ 

dr* 



■Xo 



(2r h G-r 2 h G 7r )f + r 2 h f, v + G,e 



_d_ 

dr* 



Xi , 



(17) 



in which the coefficient A is an infinite form of 0/0- type. Its value can be computed by 
means of L' Hospital rule as 



.4 



lim 

r->r h (v ,9 ) 



2r 2 (G - r h ) + r 4 / 2 + 2fr 2 r' h + r' h 
r -r h 



2r 2 h G >r + Ar h {G - r h ) + Ar\f 2 + 4fr h r' h 
2r 2 h G, r + 2r 3 J 2 - 2r' h 2 /r h . 



Now we adjust the parameter k in Eqs. (|16|]PT|) such that the ratio of the coefficients 
o ■> and 75 — o — 

or; or t ov, 

A 



before A and a d \ is equal to 1 : 2, namely 



— + 2rl(2G - f fc ) + 2rJ/' + = , 

and get the surface gravity of the event horizon 

r hG,r + r 3 J 2 - r' h 2 /r h 



K 



rl(l-2G)-riP + r' h 
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Dividing both sides of Eqs. (|16|,|l7D by r\ after substituting the coefficient A into them, 
we can reduce both of them to a standard wave equation near the event horizon 

d 2 d 2 d 2 d 

— * + 2-—-* - 25—-* + 2C— V = , (19) 

where B = f + r^/r^, for * = xo we have 

2G r' r" 

2C = — + G, r + 2r h f - (2fr h + -|) cot 9 - -\ , 

while for $ = xi we have 

2G r' 
2C = G r -(2/ + 4)cot^ + 

+ t±^!A [ {2 r h G-r 2 h G, r )f + 

From the above prescriptions, we can see that it is with the aid of the relations 
derived in the preceding section that we can recast the coupled equation (|17|) into a second- 
order wave equation of one component. This demonstrates that our disposition of the first- 
order and second-order equations is self-consistent and meets with the needs on physical 
grounds. 



2r' h 2 


r" 

— _A 






r hf,v 


+ G >6 



5 Thermal radiation spectrum 



Now that all coefficients in Eq. (fl9|) can be regarded as finite real constants, one can 
separate variables as 

* = R(rjQ(0*)e i< - mif '- uv * ) , (20) 
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and with this ansatz substituted into equation fll^D, one gets the following equation 

6' = A6, R" = 2(iu - C q )R! , 

with their solutions being 
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,A0, 



R = C ie 2(^-C )r» + ^ 



:2D 



where A is a real constant introduced in the separation variables, Co = C — XB. 
The in-going wave and the out-going wave are, respectively 



i(mifi— u!v t )+A9* 



^out ~ e i(m<p-uv*)+\$* ^(iu>-Co)r* 



:22i 



Near the event horizon, we have r* ~ =- ln(r — r\). The in-going wave \l/i n is regular at 
the event horizon, but the out-going wave \& ut(^ > ^h) is irregular at the event horizon 
r = rh, it can be analytically extended from the outside of the hole into the inside of the 
hole through the lower complex r-plane 



r - r h 



[r h -r)e 



to 



out 



(23) 



According to the method suggested by D amour- Ruffini-Sannan, 03 the relative scat- 
tering probability of the outgoing wave at the horizon is 

2 



out 



out 



-2ttlu/k 



(24) 



and the thermal radiation Fermionic spectrum of Weyl neutrinos from the event horizon 
of the hole is given by 

Wo))- -0^. (*) 
with the obvious expression of the Hawking temperature being 



Mr h — rla cos 9 



j 2 



(26) 



2tt Anr h Mr h + r\a cos 9 + r' h 2 /2 ' 

It follows that the temperature depends not only on the time, but also on the angle 9. The 
temperature is in accord with the corresponding result presented in Ref. [BJ. 



7 



6 Conclusions 



Equations (|T0| ) and ( p6|) give the location and the temperature of event horizon, which 
depend not only on the advanced time v but also on the polar angle 9. They are just the 
same results as that obtained in the discussion of the thermal radiation of scalar fields. § 
Equation fl2"o| ) shows the Hawking radiation spectrum of Weyl neutrinos in a rectilinearly 
non-uniformly accelerating Kinnersley black hole. 

This study and those in Refs. [4,12-14] manifest that the method of generalized coor- 
dinate transformation is a powerful tool to investigate black hole radiation. It was initially 
proposed by D amour- RuffiniGt to deal with the quantum thermal effect of scalar fields in 
a static (Schwarzshild) black hole and in a stationary (Kerr) black hole. Later SannanQ 
developed this method and made it suitable not only to dealing with the thermal radiation 
of scalar fields but also to that of fermions. Zhao Zheng0 further extended this trans- 
formation in the static and stationary cases to that in the non-static and non-stationary 
cases and have made several discussions about the quantum thermal effects of scalar fields 
in various black holes and that of Dirac particles in some non-static and spherically sym- 
metric black holes. But this method meets great difficulties when it is applied to the 
Hawking radiation of Dirac particles in a most general space-time such as a non-static and 
non-spherically symmetric black hole or a non-stationary axisymmetry black hole. The 
reason is that the Dirac equation can not be completely separable for all variables in the 
most general space-time. Now this difficulty is overcome when we simultaneously treat 
the first-order and second-order equations, and use the relations between the first-order 
derivatives to remove the crossing-terms in the second-order equations. This procedure 
enables us to reduce each second-order equation of one component to a simple separable 
standard wave equation near the event horizon after we have made a generalized tortoise 
coordinate transformation. 

In conclusion, the work in this paper and those in Refs. [4,12-14] make the method of 
generalized tortoise transformation more perfect. They indicate that this theory becomes a 
fairly integrated system. The method developed by us is applicable not only to discussing 
the Hawking radiation of Weyl neutrinos in a non-static and non-spherically symmetric 
Kinnersley black hole with arbitrarily acceleration and that of electrons in a non-stationary 
axisymmetry black hole, but also, in principle, to studying the quantum thermal effect of 
any space-time with a non-degenerate event horizon. Also it is easily extended to discuss 



the Hawking effect of fields with an arbitrary spin.t!-^ 115 
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